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Abstract 

This paper studies the throughput performance of HARQ (hybrid automatic repeat request) protocols over 
block fading Gaussian channels. It proposes new protocols that use the available feedback bit(s) not only to 
request a retransmission, but also to inform the transmitter about the instantaneous channel quality. An explicit 
protocol construction is given for any number of retransmissions and any number of feedback bits. The novel 
protocol is shown to simultaneously realize the gains of HARQ and of power control with partial CSI (channel 
state information). Remarkable throughput improvements are shown, especially at low and moderate SNR (signal 
to noise ratio), with respect to protocols that use the feedback bits for retransmission request only. In particular, 
for the case of a single retransmission and a single feedback bit, it is shown that the repetition is not needed at 
low SNR where the throughput improvement is due to power control only. On the other hand, at high SNR, the 
repetition is useful and the performance gain comes form a combination of power control and ability of make up 
for deep fades. 
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I. Introduction 

IN today networks, error correction is achieved by a combination of FEC (forward error correction) and 
ARQ (automatic repetition request). In classical ARQ protocols, a receiver requests a retransmission 
(sends a negative acknowledgment, or NACK) when an error is detected, and a positive acknowledgment 
(ACK) otherwise. In this work we explore the performance gain achievable by using the retransmission 
request bit(s) to signal to the transmitter the decoder status and the actual channel state, albeit coarsely. 
Our goal is to simultaneously enable the performance gain due to HARQ (hybrid automatic repeat request), 
i.e., a combination of ARQ and FEC error control methods [j3|, and to power control at the transmitter [|4|. 

A. A Motivating Example 

Consider a fixed rate transmission scheme over a block fading Gaussian channel with unit noise power 
spectral density. Let the transmit power in slot t, t E N, be Qt, the fading power gain be 7^, and the 
transmission rate be R. The receiver fails to decode when the instantaneous channel capacity is below the 
transmission rate [5], in which case it feeds back a NACK to the transmitter. A NACK is thus equivalent 
to 

- 1 

log(l + -ftQt) <R^-ft< —7^, 

vt 

that is, a NACK is a 1-bit quantization of the channel state information (CSI) 7f sent to the transmitter 
when the transmitter no longer needs it (as already remarked in [[6|, [|7)). This simple observation raises 
the question investigated in this work: whether it is optimal, in some sense, to feedback a ACK/NACK at 
the end of a slot, or whether the same feedback resources should rather be used at the beginning of the 
slot to inform the transmitter about the instantaneous channel quality, albeit coarsely. 
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To further gain insights into the problem, consider the outage capacity [j5| as the performance measure. 
For a fixed positive parameter s, let the transmission rate be parametrized as R = log(l + P s), where 
P denotes the average SNR P = E[Qt]. As explained before, a 1-bit feedback used for ACK/NACK at 
end of slot t indicates to the transmitter that 7^ < s if a NACK is received, or that 7^ > s if a ACK 
is received. The probability of successful decoding is then the probability of receiving a ACK; thus the 
outage capacity, or long term average successfully decoded rate, is 



On the other hand, consider the case where the 1-bit of feedback is used at the beginning of the slot 
to indicate to the transmitter which of the events, {7^ < s} or {7^ > s}, has occurred. In this case the 
transmitter can use this information as follows. It turns transmission off (i.e., Qt = 0) if the channel is bad 
(i.e., 7t < s) and it sends with power Qt = if the channel is good (i.e, 7^ > s); the chosen transmit 
power Qt = is such that no outage occurs. The average transmit power of this simple power control 
policy based on 1-bit CSI is P = Pr[7t > s] and the outage capacity is 



It is immediate to see that, for the same set of parameters P and s, and with 1-bit of feedback in both 
scenarios, the outage capacity with CSI ?7csi is larger than the outage capacity with ACK/NACK t^ack- 
This observation reinforces the idea that using the 1-bit feedback to signal ACK/NACK is not optimal in 
general. The question whether this conclusion changes if retransmissions are allowed is investigated in 
this paper. 

B. Past Work 

To the best of the author's knowledge, past work available in the literature considering quantized and/or 
noisy CSI only focused on outage capacity, or on outage probability, or on expected capacity, but not on 
HARQ protocols. 

For example, in [[8| the authors consider power control policies for minimizing the outage probability 
with partial CSI; the derived power policy shows benefits with respect to the case of complete absence 
of CSI even if the channel knowledge is noisy and/or partial; the benefits are more pronounced at low 
SNR. In J9|, the authors considered the outage capacity with the so called "broadcast approach", that is, 
a multiple-layer coding scheme with infinite many layers, where the receiver decodes as many layers as 
possible given the actual channel fading; it is found that even 1-bit of CSI helps to improve performance. 

In [10], the authors studied the ergodic capacity of channels with states where the state is only partially 
known at the transmitter; for the Gaussian channel with quantized CSI, they showed that the capacity 
achieving power allocation is of the waterfiUing type. In [[6|, the authors studied the expected capacity 
with quantized CSI and multiple-layer coding schemes; they showed that multiple-layer transmission offers 
limited benefits when power control at the transmitter is possible. 

In [7 1, the authors considered the DMT (diversity multiplexing tradeoff) of multi-antenna channels with 
HARQ; in this setting the feedback is only used to signal ACK/NACK and not to perform power control, 
even though the transmitter is allowed to vary the transmit power across retransmissions; it is found that 
HARQ improves the DMT by a factor proportional to the maximum number of repetitions. 

C. Contributions 

In this work we consider the joint design of HARQ protocols and power control for block fading 
Gaussian channels. We use the long-term average decoded rate Q, or simply throughput for brevity in 
the following, as a measure of performance. The throughput captures the fundamental performance limits 



Vack = Pr[7t > s] log(l + Ps). 
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when strict delay constraints are imposed and includes the outage capacity and the ergodic capacity as 
special cases. 

When considering HARQ protocols, it is customary to assume that the transmitter has no knowledge 
of the instantaneous fading and it thus transmits with equal power in every slot [|3jl. However, especially 
with INR (incremental redundancy), the probability of having to transmit m channel packets per data 
packet is decreasing with m [7]. Thus, it is conceivable that using more power in earlier transmissions of 
the same data packet reduces the probability of decoding failure and hence increases the throughput for 
the same average transmit power. Moreover, if the fading is known at the transmitter, more power can be 
used in the most favorable channel conditions-assuming power control is possible. As pointed out in [jTT|, 
in delay constrained scenarios, the assumptions about the dynamics of the fading process with respect to 
the code length, as well as the duration over which power constraints are enforced, are critical. Here, in 
order to enable power allocation, we consider a power constraint imposed over a time horizon comprising 
many slots (i.e., much bigger than the maximum number of retransmissions allowed), commonly referred 



to as long-term average power constraint [12|. 



As opposed to classical HARQ protocols, where the ACK/NACK feedback bit is sent at the end of the 
slot, we consider here systems where the feedback bit(s) can be sent back at any point in time during 
a slot; the feedback can be used to signal CSI, or ACK/NACK, or any combination of them. The only 
restriction we impose is that the feedback bit(s) cannot be carried over from slot to slot. 

Our contributions can be summarized as follows: 

1) We propose novel HARQ protocols where the CSI and the ACK/NACK information are combined 
within the same feedback bit(s) in order to realize simultaneously the gains due to HARQ and the 
gains of power control with partial CSI. The main idea behind the proposed protocols is that the 
receiver sends back to the transmitter the index of the smallest power level that will allow successful 
decoding in the current slot. Our protocols are time-varying quantizers for a suitably scaled version 
of the channel fading, where the scaling factor accounts for the information already available at the 
receiver from the past transmissions. 

2) We show that the throughput performance of the proposed class of protocols with perfect CSI can 



be obtained from dynamic programing [13| 



3) By numerical evaluations of the throughput for Rayleigh fading channels, we show that repetitions 
are not needed at low SNR, and that the improvement over classical HARQ protocols (that use the 
feedback bit for ACK/NACK only) is entirely due to ability to perform power control. 

At high SNR, repetitions are useful, and the performance improvement over classical HARQ proto- 
cols comes form a combination of power control and the ability to make up for atypical long deep 
fades with repetitions. 

Our numerical results show that our protocols outperform classical HARQ at all SNRs. 

4) We also have the following side results: (a) we show that the optimal power allocation for the 
outage capacity with partial CSI consists of a quantizer of the fading gain where the quantization 
regions are union of intervals, rather than intervals; to the best of the author's knowledge-was not 
reported before; (b) we present novel bounding techniques for to compute certain probabilities that 
are needed for the throughput evaluation; these techniques are useful for numerical optimizations. 
In particular, a technique based on considerations on the order statistics of an independent sample 
of negative exponential random variables is of interest in its own. 

D. Paper Organization 

The rest of paper is organized as follows. Section |ll] introduces the system model and Section III 



evaluates the throughput; Section |IV] revises the ergodic capacity with partial CSI, which serves as an 
upper bound for any HARQ protocol; Section [V] derives the outage capacity with partial CSI, which serves 



as a lower bound for any HARQ protocol; Section VI proposes a new class of HARQ protocols that 



combine repetitions and power control for any number of retransmissions and any number of feedback 
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Fig. 1. The block-fading Gaussian channel model with partial CSI at the encoder (dotted blue line) and perfect CSI at the decoder (solid 
red line). 



bits; Section |VII| proposes a novel bounding technique for the throughput based on ideas from order 
statistics; Section |VIII compares the throughput performance of our new protocols with that of classical 
HARQ protocols for the Rayleigh fading channel; Section |K] concludes the paper and points out open 
questions and future work directions. 



II. System Model and Performance Metric 

We adopt the following notation convention: [x]"*" indicates max{x, 0}, l{xeA} the indicator function 
(that equals one whenever x E A and zero otherwise), and Fx{x) = Pr[X < x], x G M, is the cumulative 
distribution function of the random variable X. In the following slot, fading block and codeword length 
are used interchangeably. 

This work considers the single-user block-fading Gaussian channel. The received signal vector in slot 

t,teN, is 

Yt = ^/iTQtXt + Zte C^ 

where: the noise Zt is a length-L proper-complex white Gaussian random vector with zero mean and unit 
variance, the channel fading power gain 7t is a scalar with E[7t] = 1, the channel input signal Xt has 
Gaussian iid (independent and identically distributed) components with zero mean and unit variance, and 
Qt >0 is the transmit power. Each codeword Xt spans one fading block over which the fading gain stays 
constant. The slot length L is such that it suffices to guarantee reliable communication if the accumulated 
mutual information at the receiver is above the communication rate [j5|. 

The fading gain 7^ changes in an iid fashion from slot to slot. The receiver has perfect instantaneous 
knowledge of 'jt at the beginning of the slot. The transmitter however does not know 74, unless explicitly 
informed by the receiver. For this reason, we assume that the transmitter cannot adjust the communication 
rate in each slot and thus it sends at a fixed rate. The (partial) CSI possibly available at the transmitter is 
used only for power allocation across slots 

A delay-free and error-free feedback channel with capacity log2(-F) bits per slot is available for 
communication of low-rate information between the receiver and the transmitter; the receiver can feedback 
a retransmission request to the transmitter at the end of a slot, or quantized CSI at the beginning of a 
slot, or any other information representable on log2(-F) bits at any point during the slot. We do not allow 
feedback bits to be accumulated over successive slots. The case F = 1 corresponds to absence of CSI at 
the transmitter, while F = +00 corresponds to perfect CSI. 

In a block-fading setting, reliable communication is possible if the accumulated mutual information at 
the receiver is above the communication rate [|5J. To make up for decoding errors, which occur when 

'We note that in practice, fading can be considered independent from slot to slot only if the slots are separated in time by at least few 
channel coherence times [14J . This is not a problem in multi-user systems where a user is assigned a transmission slot in every frame (and 
frames consist of several slots). In this case 72 indicates the fading gain of the slot where the second transmission (first re-transmission) 
occurs and it needs not be the slot immediately following the one where the first transmission occurred. When the iid assumption does not 
hold, fading correlation across time slots can be easily incorporated in our model by substituting products of probabilities involving different 
fading random variables with the corresponding joint probabilities. 
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the channel is in deep fade, the transmitter can retransmit a data packet at most M — 1 times, that is, 
each data packet can be transmitted on at most M channel slots. We consider the three HARQ protocols 
analyzed in [|3J: 

• ALO (ALOha): like in slotted Aloha, the transmitter keeps sending the same codeword and the 
receiver attempts decoding by using only the most recently received codeword. 

• RTD (Repetition Time Diversity): the transmitter keeps sending the same codeword and the receiver 
performs maximal ratio combining of all the received packets, thus realizing Repetition Time Diver- 
sity. 

• INR (INcremental Redundancy): at each retransmission request, the transmitter sends new redundancy 
bits and the receiver optimally combines them. 

The protocols work as follows. In order to send a data packet of b bits, the transmitter can use at most 
LM channel uses. For ALO and RTD, the transmitter encodes the data packet with a Gaussian channel 
code of rate b/L and then concatenates it with a repetition code of rate 1/M. For INR, the transmitter 
encodes the data with a Gaussian channel code of rate b/{LM) and at each transmission it send a different 
chuck of L symbols. Following [jsj, we let i? = The throughput is defined as the long-term average 
number of successfully decoded bits per channel use: 

1 ^ 

^M,F = — y R ■ 1 {successful decoding in slot t}j (1) 



t= 



where the subscript M indicates that at most M transmissions are allowed per data packet, and the 
subscript F indicates that at most F feedback values are allowed. The ergodic capacity and the outage 
capacity are a spacial case of our framework for M = +00 and M = 1, respectively. For future use, 
decoding fails with m transmissions (and thus a retransmission is needed) if: 

m \ 

J2 log(l + ItQt) <r \ for INR (2a) 
t=i J 

log(l + ItQt) <r \ for RTD (2b) 
t=i J 

U I log(l + ^tQt) < r] for ALO, (2c) 



t=i 



since INR accumulates mutual information, RTD accumulates SNR, and ALO only accounts for the most 
recent transmission. 

In order to complete the system description, we need to specify how the transmit power Qt, t G N, 
can be varied. We assume Qt G {Pmj}, where Pmj > is the power used at the m-th transmission 
attempt, m G {1, . . . , M}, when the feedback value is /, / G {0, . . . , F — 1}. The total transmit power 



must satisfy the long-term average constraint / [72] / defined as: 

T 



lim ^ T^Qt < P-, almost surely. (3) 



t=l 



For the normalization adopted in this work, P has the meaning of average SNR at the receiver. The 
power allocation policies are causal pT| , p3| in that Qt can only depend on the (partial) knowledge of 
(7i,---,7i)- 

Next we show how to evaluate the throughput in ^ subject to the power constraint in Q. 
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III. Throughput Evaluation 
In [|3| we introduced a general framework to analyze the performance of HARQ protocols based on the 



renewal-reward theory [16|, which we shall use now to evaluate the throughput. In our system, a renewal 
event occurs (i.e., the system starts anew without any memory from the past) when the transmission of a 
data packet ends (either because of successful decoding with less than M transmissions or because the 
maximum number of transmissions has been reached). From ||3|, the system performance is completely 
characterized by the triplet (T,TZ,V), where: 

• T G {1, . . . , M} in the inter-renewal time and represents the number of slots needed to complete 
the transmit of a data packet; 

• 7^ G {0, R} is the reward, i.e., the number of bits successfully decoded per channel use when the 
transmission of a data packet ends; 

• P is the cost, i.e., the total transmit power for a data packet (including all retransmissions); 
Given {T,TZ,V), the throughput in ([T]) subject to the power constraint in Q is given by: 

Theorem 1 (from [[3|, [|7|). For any M and F the throughput rjM f far a given power P is the solution 
of: 



Vm,f 



E[n] R{1-Ft[T = M, failure to decode]) 



max —T—r = max 



^ _ Eili nVmir > m] Pr[r > m] ^ - 

"nn E™=iPr[r>m] - ' 

where the maximization is over the transmit rate R> and over the power allocation {Vm}m=v where 
Vm G {-Rm,0! • • • 5 -fm,F-i} the causal power policy far the m-th transmission attempt restricted to take 
at most F different values. The distribution of the inter-renewal time T and the probability of failure to 
decode (see (|2])j are function of the protocol ^ G {ALO, RTD, INR} used. 

Proof: The proof can be found in Appendix |Aj ■ 
Remarks: 

1) The probability of failure to decode on the last transmission, indicated as Pr[T = M, failure to decode], 
is the probability that the data packet is lost and it is referred to as outage probability. 

2) The event {T < M — 1} (that transmission ends before the maximum number of transmissions has 
been reached) implies successful decoding. However, a successful decoding does not necessarily 
imply a renewal event (the end of the transmission of the current data packet). 

3) It is immediate to see that the optimal power allocation meets the power constraint with equality 
(otherwise, the left over power could be used on the last transmission, which would increasing the 
throughput while still meeting the power constraint). 

The throughput of the different protocols for different values of M and F satisfies: 

Theorem 2 (from f3]). We have: 

.sfp°' < .ir' < .c''- (■») 

Moreover, rj^^jp is a non-decreasing function of M and of F, for each protocol -k G {ALO, RTD, INR}. 

Proof: The proof of Q is as in (3\ and is omitted here for sake of space. The fact that rj'lfp is a 
non-decreasing function of M and F follows by observing that if more transmissions or more accurate 
CSI would hurt performance, they could just be ignored. ■ 
From Theorem [2] it follows immediately that: 
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Corollary 3. For any M and F, and for any protocol * e {ALO, RTD, INR}: 

(ALO) ^ {*) ^ (INR) 
Vm=1,F — Vm,F — %/=oo,F' 

where rf^^^'p is the outage capacity of the channel and rf^^^ p is the ergodic capacity of the channel 

In the following we first evaluate ^7^1^ p and with partial CSI and then we propose novel 

achievable protocols for r/^^^^,, e {ALO, RTD, INR}. 

IV. Throughput upper bound r/^™^^^ 

When M — !■ oo, the INR protocol with a time-invariant and memoryless power allocation policy 
Pt = gilt), t E N, and with optimized rate R, achieves the ergodic capacity of the channel [3| given by: 

CJE,F = E[log(l + 7^?(7))], 

where the function g{-) can take at most F different values. In the following we let Pj > be the 
power used when the feedback value is /, / G {0, . . . , F — 1}, (i.e., we drop the index referring to 
the number of repetitions, which is irrelevant here because we considered time-invariant and memoryless 
power allocation policies). The optimal power policy g(-) for a finite F was derived in [lOj and it is 
summarized in the following: 

Theorem 4 (from [[TOj). Let 

0<Po<Pi<---<Pf-i<j, (5) 
where Pf is the power used when 7 G 7l^j^^\ with 

7^P = W<7<^/+l}, (6) 



/ G {0, . . . , F — 1}, with So = 0, sp = +00, and 



1 A 1 ( - Pf) .p , .7. 



/ G {1, . . . , F — 1}, and where X> is such that 

Xip/Pr[7^P^] = P. (8) 

/=o 

The ergodic capacity with partial CSI is: 

(INR) _ 
Vm=oo,F — 
F-l 

= max^E[log(l + ^Pf)\-f G TZf^] Pr[7^5.™^^], (9) 

where the maximization is subject to ([5]) and ([8]). 
Remarks: 

1) When F = 1 (no CSI) the transmitter can not adapt its power across transmissions and thus sends 
at constant power (7(7) = P. In this case the throughput is: 

r/E,^=i = E[log(l + 7P)]. (10) 
Theorem |4] for F = 1 gives the result in ( flQ] ). 
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2) With F = oo (perfect CSI) the optimal power allocation is water-filling [|4) given by: 

1 1' 



and the throughput is: 



(INR) 
Vm=oo,F=oo 



X 7 



E 



log 



7 
AJ 



(11) 



(12) 



where the Lagrange multiplier A > is such that the power constraint is met with equality. The 
power policy in (|7]) for F — )• 00 reduces to ( [TT] ) since ^pJ'^^^Sp^)^ ^ — )■ 1 when Pf^i — Pf — t- (notice 

that the region TZq^^^ always includes the interval [0, A] since Si > A, while the other quantization 
regions reduce to a single point). 

3) The optimal quantization regions are intervals. 

4) For the purpose of numerical evaluations, it is convenient to have bounds on the throughput that 
can be fast evaluated and easily optimized. The throughput in Theorem |4] can be bounded as: 
Proposition 5. For a given set of quantization intervals IZf = {s/ < 7 < / G {0, . . . , F — l}, 

(INR) 



the ergodic capacity tim^Jo f ® bounded as: 



F-l 



^r=^i,F < max E log (1 + f^f Pf) Pr[7^;], 



7=0 

F-l 



C=^i,F > max V log(l + SfPf) Pr[7^y], 



(13a) 



(13b) 



where fij = ]E[7|7 G TZj] 
maximization is subject to 



G 71 f is the centroid of the f-th quantization interval, and the 



= So < si < S2 < ■ ■ ■ < sf = +00 



(14) 



and such that the powers Pf >0, / G {0, . . . , F — 1}, satisfy X]/=o 'P^l ^ T^f] ^ P- 



Proof: The bounds in ( [T3| ) follow immediately from the definition of the quantization intervals 
and from Jensen's inequality, i.e., 

log(l + inf{7G7^/}P/) 
<E[log(l + 7P/)|7G7^/] 
<log(l + E[7|7G7^/]P^), 



with inf{7 G 7^/} = s/ by definition. For both bounds in ( [T3] ) the optimal powers are obtained by 
water-filling [4|. The optimization of the bounds in ([T3]) is thus equivalent to the problem of finding 



the optimal quantization intervals, which can be done efficiently by using Lloyd's algorithm [17| 



As the number of feedback levels F increases, the quantization intervals reduce to a single point 
since fxj — )■ sj, that is, the bounds in ( [T3] ) converge to the water- filling ergodic capacity in ([12]). 
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V. Throughput Lower Bound r7it=?,F 
All protocols have the same throughput for M = 1 (because retransmissions are not possible), which 
coincides with the outage capacity of the channel [|3| given by: 

Vuti^F = max {R{1- Pout(i?))}, 
where Pout(-R) is the probability of outage given by: 

Pout(i?) = Pr[log(l + 7^7(7))<i?], 

and where the non-negative function g(-) can take at most F different values. The optimal power allocation 
policy g(-) for a finite F is: 

Theorem 6. Define the thresholds 

< si < ■ ■ ■ < sp-i < sp = So < sp+i = +00, (15) 

(notice the convention sp = Sq & (0, oo), rather than sq = and sp = +oo as in ([T4])j and the 
quantization regions: 

^(ALO) = {0 < 7 < U {7 > So}, (16a) 

)(A 



7^(ALO) ^ < ^ < / G {1, . . . , F - 1}, (16b) 



Let the transmit power to be used when 7 G TZ^^^'^^ be: 



Pf = ^!^, /G{0,...,F-1}, (17) 

■5/ 



r/ze outage capacity with partial CSI is: 



'IM 



= max ^1 — Pr[7 < Si]j- 



where the maximization is subject to the constraint in ( [15] ). 

Proof: The proof can be found in Appendix |B} 
The power Pf in ( [17] ) is the minimum power that guarantees no outage for all fading gains in 7^^^^*^^ 
for / > 0; outage can only occur when the fading gain belongs to the subset of T^-q^^*^^ given by {7 < si} 



Remarks: 

1) When F = 1 (no CSI) the transmitter can only send at constant power (7(7) = P and the throughput 
is: 

vlit?F=i = maxlog(l + Psi) Pr[7 > s,]. (19) 

' si>0 



Theorem [6] for F = 1 gives the result in ( [19] ). 



2) With F = 00 (perfect CSI) the optimal power allocation is truncated channel inversion [ [T2| : 

- 1 

^(7) = — — l{7>si}- (20) 



10 



With (|20l) an outage only happens when 7 < si, and the throughput is 



(ALO) 
Vm=1,F=oo 



maxPr[7 > Si] log 

si>0 



1 + 



E 



(21) 



7 -L{7>si} 



When F — )■ 00, the region TZq 

,(ALO) 



(ALO) 



reduces to {0 < 7 < si} since Soo = sq = 00; when the 
fading belongs to IZ)^"^^' the transmit power is Pq = Vmis-^oo ^-f^ = 0, thus Theorem 6 reduces to 
truncated channel inversion. 

3) Several power allocation policies have been proposed for outage minimization with partial CSI. 
However, none of the policies is optimal. For example the solution proposed in [8] corresponds 
to the suboptimal solution sq = sp = sp+i = 00, that is, setting the quantizations regions to be 
intervals. From our result in Theorem [6} the optimal quantization regions are in general unions of 
intervals. 

4) For the purpose of simplifying our numerical evaluations we propose to bound the throughput 

(ALO) 
11m=1,F 



Proposition 7. Let 



;^(ALO) 

'Im=i,f 



■ log 1 



max Pr[7 > Si] 

0< si < . . . < Si? < s j7_|_i =00 

P 



(22) 



The throughput in (18) is bounded by 



ij(ALO) 

'/m=i,f-i 



< n[f°' 



< ^ALO) 
— 'IM 



l,F- 



When F 1, the optimal solution of ( [22] ) tends to 

Sf = s,e~\Wf>l, 

for some ^ > 1. 



(23) 



Proof: That rfjf^ll^^p in ( [22| ) is an upper bound for ri)^^i'p in ( (T8] ) follows by neglecting the 
term ^''^J^''^^ at the denominator in ( [T8] ). That rfj^^^^p_^ in ([22]) is an lower bound for r^^'^^f ^ in ( [T8| ) 
follows by setting Si? = 00 in ( [T8[ ). It is interesting to notice that the same function r^^^f^ is a 



(ALO) 



lower bound for r7^=?F+i (^^otice the different number of feedback values) and an upper bound for 

„(ALO) 

'/m=i,f- 

The proof of ( 123] ) can be found in Appendix [C] ■ 



VL Main result: achievable throughput for general M and F 

Determining rj'i^p, G {ALO, RTD, INR}, as in Theorem |l] for general finite values of M and F, is 
very complex as it involves the solution of a dynamic program (due to the causal nature of the power 
control pT| , [15|). In this section we propose novel protocols that combine repetition and power control 



for a general M and F. The throughput of our protocols is a lower bound for the optimal tj^mf- 

Theorem 8. For a protocol -k G {ALO, RTD, INR} and general finite values of M and F, let Bm G 
{0, . . . , -F — 1} be the feedback sent by the receiver at the beginning of the slot corresponding to the m-th 
transmission attempt for the current data packet, m G {1, . . . , M}. Consider the following power policy: 
form G {1,...,M} 



Trr 



-{B„ 



F-1 

=o} + E 

/=i 



- 1 



1 



{Bm=/}' 



(24) 
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with 



< r„ 

= Smfl < Sm,l ■ ■ ■ < Srri,F-l < Sm,F = +00. 



(25a) 
(25b) 



The thresholds {smj}f=o in ( |25| ) define a quantizer for a scaled version of the fading power gain 7^, where 
the scaling factor accounts for the information already accumulated at the receiver in the previous m — 1 
transmissions. In particular, the proposed feedback policy is: for -k G {ALO, RTD, INR}, m G {1, . . . , M} 
and f e {0, . . . ,F - 1} let 



Bm = f if ^ e [smj, Smj+i) and ^l*> > 0, 



Br, 



for with = 1 and for m > 1 defined as 



(ALO) 
m 



^•(RTD 
(r(INR) 



m—1 

ni{i-^>o}^ 

m—1 



E 



t=i 



nr=7(i + (e^-i)S 



(26a) 
(26b) 

(27a) 
(27b) 

(27c) 



The resulting throughput for protocol -k G {ALO, RTD, INR}, as given in Theorem^ is lower bounded 
by: 



i-k) 1 ~ Pout 

VM,F,lh — . , T , v^M-1 



^m=l 
^ + Z^m=l Pni,0 



■ I02: 1 + P 



(28) 



m=l 2^f = l Sm,f 

where the maximization is subject to the constraints in ^25^ and where the probabilities {pmj} ^re defined 



as: 



-{*) 

PmJ = PlnJ 



m = l,...,M, / = 0,...,F-1, 
{by defining sm+i,i = +00), 

where {p^ are defined in ( [361 ) fat" ALO, in ([37]) for RTD, and in ( [38] ) for INR. 



P - 

r out — FM+1,0 



(29a) 
(29b) 



The rest of the section is devoted to give a rational for the protocols in (|26])-([27]), to prove the throughput 
formula in ([28]) and to define the probabilities {p^y} for ([29]). 



A. Protocol description and throughput evaluation 

Inspired by the power policy that minimizes the outage capacity in Theorem [6] we propose that the 
transmitter uses the power policy in ( [24| ). Our protocol works as follows: 

• The receiver feeds back Bm = / > 0, m G {1, . . . , M}, to indicate that the power (e^ — l)/smj 

suffices to successfully decode the current data packet when the previous m — 1 transmissions are 

combined with the current transmission. 
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• Upon receiving Bm = / > 0, m G {1,...,M}, the transmitter is certain that the receiver will 
decode correctly with the current transmission; hence, after transmission with power (e^ — 

the transmitter prepares to sent a new data packet. 

• The receiver feeds back B^ = 0, m G {1, . . . , M}, when none of the powers (e^ — l)/smj, V/ > 0, 
would guarantee successful decoding. 

• In response to Bm = 0, m G {1, . . . , M — 1}, the transmitter sends with power (e^ — l)/rm and 
prepares to retransmit the same data packet in the next slot. 

• In response to Bm = (for the last transmission attempt), the transmitter sends with power (e^ — 
1)/tm and prepares to sent a new data packet since no more retransmissions are permitted. In this 
case an outage can occur. 

In order to evaluate the throughput according to Theorem [1] we must determine the average decoded 
rate (reward) and the average transmit power (cost) when the transmission of the current data packet ends, 
and the average time needed to transmit a data packet (inter-renewal time). From the description of the 
protocol given above, it is clear that the transmission of the current data packet does not end after m 
received slots, m G {1, . . . , M — 1}, if all the feedback values received were zero, that is, 

Pr[r >m] = Fi[B, = 0, . . . , = 0]. 

When the transmission of a data packet ends with less than M transmissions, successful decoding occurs. 
The transmission of the current data packet ends after the M-th transmission regardless of the status of 
the decoder (since no more transmission attempts are possible). An outage occurs if decoding is still 
unsuccessful with M transmissions, i.e.. 

Pout = Pr[5i = 0, . . . , Bm = 0, failure to decode]. (30) 

The average number of successfully decoded bits when the transmission of the current data packet ends, 
i.e., average reward, is 

E[7^] =i?(l-Po,t). (31) 
The average transmit power when transmission of the current data packet ends, i.e., average cost, is 

^[^] = P^t^i = 0, . . . , 5^ = 0] 

m=l 



M F-l 



+ 1^1^ P^t^i = 0' • • • = 0'^™ = /]■ (32) 



m=l /=1 ^"^'f 



The average time needed to transmit a data packet, i.e., the inter-renewal time, is 

M M 

E[71 = Pr[r = m] = ^ Pr[r > m] 



m=l m=l 
A/-1 



l+^Pr[5i = 0,...,5^ = 0]. (33) 



m=l 



The probabilities in ( [32| ) and (|33j) can be easily expressed as a function of pmj defined as: 

Pmj = Pr[5i = 0, . . . , Bm-i = 0, Bm < /], (34) 

for / G {0,...,F- 1} since: 

Pr[Si = . . . = Bm-l = 0,Bm = f] 
Prn,0 ^ / = 0, 

PrnJ - PmJ-l / = 1, . • . , F - 2, (35) 

Pm-1,0 / = F - 1, 
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for all m G {1, . . . , M}. The equality for / = F — 1 follows since Pr[_Bm < — 1] = 1 for all m. As we 



shall proof in the next sections, where we give the details of the protocols, the outage probability in pO| ) 

_i if we assume an hypothetical 1-bit of feedback at the end 



l{failure to decode with M transmissions}- (this wiU COrrCSpOnd tO a 



needed for pT] ) ie equivalent to Pout = Pm 
of the M-th transmission that indicates B^+i 
degenerate quantizer with = sm+ia < sa/+i,i = +oo). This discussion justifies the definitions in ([29 
Finally, with the definition in ( [35] ) and by Theorem [T] the throughput for the proposed protocols is 

_ eq.([3T]) 



eq.(33) 



eq.(28) 



where we expressed the rate R in ( [31] ) as a function of P from 



eq.( 


32) 


eq.( 


33) 



-1)- 



Pm.O _|_ Pm,f-Pm,f~ 



1 + E 



M 

m=l 



Pmfl 



The probabilities {p^ in ( |29[ ), as well as the feedback policy feedback policy {Bm} in ([26]) defined 

through the quantities {Qn^ } in ([27]), depend on the protocol -k G {ALO, RTD, INR} and will be discussed 
next. 



B. First transmission (all protocols) 

In order to better understand the way the feedback value is decided, consider the transmission on the 
first slot, which is the same for all protocols. Let 9 = — 1 > 0. 
• Feedback policy: 

- At the beginning of the first slot/transmission, the receiver measures 71 and sends 

El = F - 1 if log (1 + 9^^] > log(l + 9) 

that is, the receiver sends back the highest possible feedback value if the lowest possible power 
9/si^F-i (from ([24]) with m = 1 and f = F — 1) suffices for successful decoding given the 
actual fading 71. In other words, the receiver sends back 

Bi = F — 1 if Si^F^i < 7i. 

- If Si,F-i > 7i5 the lowest possible power is not enough to guarantee successful decoding, the 
receiver checks whether the second lowest available power 9 / suffices for correct decoding. 
The receiver sends back 

5i = F - 2 if log (1 + ^^^1 < log(l + 9) 

and log (1 + I > log(l + 9), 

that is, 

5i = F - 2 if < 7i < Si.F-i- 

- By continuing our reasoning in this manner, the receiver sends 

5i = / if 7i e [si,/, / = 0, . . . , F - 1, 



i.e., the thresholds {si /} J^q define a quantizer for 71 as in ( |26a[ ) with ^0 = 1 (there is no scaling 
for the fading value on the first transmission because there no accumulate information at the 
receiver). 
Transmission strategy: 
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72 



B, = 0, So = 1 



Bi = 1 



*2,1 



B, = 0, = OK 



Bi = 0, B2 = gfi 



Tl 



71 



Fig. 2. Feedback values for the ALO protocol with M — 2 retransmissions and F — 2 feedback values. The shaded region corresponds 
to an outage. 



- In response to _Bi = / > the transmitter sends with power 0/sij that guarantees successful 
decoding for the whole range of fading values in [sij, sij+i); after transmission, the transmitter 
prepares to send a new data packet. 

- In response to Bi = 0, the transmitter sends with power 6'/ri |^ that suffices for successful 
decoding only if 71 > min{ri, si i}; however, the transmitter can not know whether the actual 
7i is above or below min{ri, and hence prepares to retransmit the same packet again. 

From the second transmission onwards, the mode of operation depends on the protocol used. We will 
describe the three protocols separately. 



C. Retransmissions for ALO 

Recall that a second transmission is triggered by i?i = 0, which corresponds to having sent with power 
9/ti on the first slot. 
• First retransmission: 

- In the ALO protocol only the most recent received slot is used for decoding. 

Assume that in the first transmission the fading satisfied 71 > min{ri, Si 1}. The receiver knows 
that the transmitter will resend the same data packet in the second slot because it received 
-Bi = in the first slot. The receiver can "trick" the transmitter into believing that it will be 
able to decode in the second slot by sending B2 = F — 1. 

Clearly, this second transmission is a waste of power, but the receiver has no other way to inform 
the transmitter of its successfully decoding owning to having already exhausted all its feedback 
bits at the beginning of the current slot. Among all possible powers that the receiver could have 



requested for the second transmission, 6/s2,f~i (from ([24]) with m = 2 and / = F — 1) is the 
lowest. This is captured in our protocol definition by the condition in ( |26b| ). 
- If 7i < min{ri, Si 1} (which implies 71 < ri), the receiver uses the same feedback policy it had 
used as on the first slot, but with possibly different thresholds for the quantization. 
Fig. |2] shows the feedback values for the ALO protocol with M = 2 retransmissions and F = 2 
feedback values; the region with Bi = B2 = is divided into two parts, the shaded region 
corresponds to an outage while the white region corresponds to successful decoding. 
Other retransmissions: 

By continuing our reasoning as for the first retransmission, we arrive at the protocol definition in 
Proposition [8] with the fading scaling defined in ( |27a| ). 



"^The power in this case is 9/ri; notice the use of ri in place of si^; this is because si,o ~ has been already used to indicate the 
left-most quantization value. The same holds for any m > 1. 
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Performance: 

For ALO the probabilities Pmj, f ^ {0, ■ ■ ■ , F — 2}, as defined in ( [34] ) are: 

= Pr [-ft < 7t_i < n_i Vt < m, 

7m-l < Tm-1, 7m < ^mj+l] 

7t 



Pr 



min{rt,Sf,i} 



< 1 Vt < m, 7m < Smj+i 



m— 1 



F^(min{rt, St,i}) j F^{smj+i), (36) 

with nt=i(' ■ ■ ) = 1- With the definition of p^j*^^ in ([36]) we obtain the relationships in ^29\ for the 
ALO protocol. 

In general, the probability p^^^^ is available in closed form if the cumulative density function of 
the fading power F^(x), a; > 0, is known in closed form. 

Notice that the probability of outage (failure to decode at the last transmission) would be equivalent 
to sending Bm+i = 0, for this reason we have the equahty in ( |29b| ) by defining sm+i,i = +oo. This 
observation holds for all protocols. 

D. Retransmissions for RTD 

Recall that a second transmission is triggered by Bi = 0, which corresponded to having sent with 
power 9/ti on the first slot. 
• First retransmission: 

- In the RTD protocol, the receiver accumulates SNR. 

On the second transmission, at the beginning of the slot the receiver measures 72 and checks 
whether the lowest available power 9/s2,f-i suffices for decoding, and sends 

^2 = F - 1 if 

log (1 + 9^ + 9^^] > log(l + 9), 

V Tl S2,F-1 J 

since the resulting SNR after maximal ratio combining of the two received packets is the sum 
of the SNR on each packet. Hence, the feedback value at the beginning of the second slot is 

^2 = F - 1 if — + > 1. 

Tl S2,F-1 

If 7^ > 1, then B2 = F — 1. However, if ^ > 1 a retransmission was not necessary in the first 
place and the power 9/s2,f-i is wasted, as for the ALO protocol. When 71 > ri the receiver 
feeds back B2 = F — 1 so that transmission of the current data packet ends with the next slot 
(with the minimum possible amount of wasted power). 

If 7i > Tl and B2 = F — 1 then the condition on the fading value we can rewrite as: 

B2=F-2ii ^^>S2,F-l. 

Tl 



- If instead — H < 1 (which implies that — < 1) the receiver checks whether the second 



2i H 2^ < 1 (which implies that ^ 

Tl S2,F-1 ^ ^ Tl 

lowest available power 9 / S2,f~2 suffices for decoding, and sends 

^2 = F - 2 



if log ( 1 + + 9^^ ) < log(l + 

Tl S2,F~l. 

and log ( 1 + + 9^^^ ] > log(l + 

Tl S2P-2' 
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Fig. 3. Feedback values for the RTD protocol with Af = 2 retransmissions and F = 2 feedback values. The shaded region corresponds to 
an outage. 



Hence, the feedback value at the beginning of the second slot is 



Bo 



F - 2 if S2F-2 < 



72 



71 

n 



< S2P-1- 



- By proceeding with this reasoning, we see that the thresholds {s2j}^=q define a quantizer for 
when ^ < 1. The value to be quantized is larger than the actual fading 72 as it 

accounts for the SNR already "harvested" at the receiver during the first transmission. 
As an example. Fig. |3] shows the feedback values for the RTD protocol with M = 2 retrans- 
missions and F = 2 feedback values; the region with _Bi = i?2 = is divided into two parts, 
the shaded region corresponds to an outage while the white region corresponds to successful 
decoding. 
Other retransmissions: 

In general, if the m-th transmission is required, then the receiver has already accumulated an 
equivalent SNR of 



SNR' 



m—l 



m—1 

t=i 



from the previous m — l transmissions, all in response to a zero-value feedback value. If SNR' 



m— 1 



> 



1 decoding was successful, else a retransmission is needed and the receiver uses the thresholds 
{sm,/}j=o define a quantizer for -^_g|^'^, — as proposed in Proposition M (fading scaling defined 
in ^^). 
Performance: 

For RTD the probabilities Pmj, / G {0, . . . , F — 2}, defined in ( [341 ), are: 



PmJ 



Pr 



It 



1-E 



7^ 



t-1 

— > Vt < m, 7„ < Smj+i 

— Tf 

1=1 ^ 



Pr 



t-1 



£=0 



7^ 



It 



min{rt, St,i} 



<l\/t <m, 



m-l 

Ti 



< S 



m,f+l 



(37) 
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with 7o/to = 0. With the definition of p^^^'^ in ( [37] ) we obtain the relationships in ([29]) for the RTD 
protocol. 

In general, there is not a closed form expression available for p^^™\ unless it is possible to evaluate 
the density of random variables of the type X]^=i 7^ closed form. 

E. Retransmissions for INR 

Recall that a second transmission is triggered by -Bi = 0, which corresponds to having sent with power 
6/ti on the first slot. 
• First retransmission: 

- For the INR protocol, the receiver accumulates mutual information. 

At the beginning of the second slot, the receiver measures 72 and sends 

B2 = F-1 if 

log (^1 + + log (^1 + Oj^^ > log(l + 9), 

since the resulting accumulated mutual information at the receiver after optimal combining of 
the two transmissions is the sum of mutual information of each slot. 

Again, if ^ > 1 then B2 = F — 1 and a retransmission was not necessary in the first place and 

the power 9/s2,f~i is wasted, as for the ALO and RTD protocols. 

If 7i > Tl and B2 = F — \ then the condition on the fading value we can rewrite as: 

1 + ^^ 

^2 = F - 1 if 72 > S2,F-1- 

Tl 

- If log ^1 + +log ^1 + 9-^^^^ < \og{l + 9) (which implies that ^ < 1) the receiver checks 
whether the second lowest available power 9 / S2,f-2 suffices for decoding, and sends 

B2 = F-2 

if log (l + 9^\ + log (l + 9^^\ < log(l + 9) 



and log ( 1 + 9 — ] + log ( 1 + 9^^] > log(l + 

V ^1 / V S2,F~2 ) 

Hence, the feedback value at the beginning of the second slot is 

1 + ^21 



B2 = - 2 if S2,F-2 < 72 ^ < S2,F-1 



- By proceeding with this reasoning, we see that the thresholds {s2 /}/=o define a quantizer for 

72 when 7i < tx\ this scaled version of 72 accounts for the mutual information already 

T\ 

accumulated at the receiver in the first transmission. 

For example. Fig. |4] shows the feedback values for the ALO protocol with M = 2 retransmissions 
and F = 2 feedback values; the region with Bi = B2 = is divided into two parts, the shaded 
region corresponds to an outage while the white region corresponds to successful decoding. 
Other retransmissions: 

For a general m, the mutual information already accumulated from the previous slots is 

m—l 



L-i = 5^1og 1 + 0^ 

t=l ^ * ^ 
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Fig. 4. Feedback values for the INR protocol with M = 2 retransmissions and F = 2 feedbaclc values. The shaded region corresponds to 
an outage. 



If 'm-i > log(l + 0) then decoding was successful and the receiver sends = -F — 1 to end 
transmission in slot m. If < log(l + 9), then the receiver looks for the smallest / > such that 



\m-l + log 1 + ^ 



7« 



> log(l + 



If such an / > exists, then B„i = f and transmission ends with the current slot; otherwise, B„i = 
and transmission continues. This procedure is equivalent to the protocol in Proposition [8] with the 
fading scaling defined in ( |27c[ ). 
• Performance: 



For INR the probabilities Pmj, / ^ {0, . . . , -F — 2}, defined in ([34]), are: 



Pr 



f(INR) 



7t 



min{Tt,st^i} 



< 1, 



(38) 



with 7o/to = 0. With this definition of p^^^'' we obtain the relationships in (|29]) for the INR protocol. 
In general, there is not a closed form expression available for P)^^\ unless it is possible to evaluate 
the density of random variables of the type 11^=1 (^1 + closed form. 



F. Performance with perfect CSI at the transmitter 

In order to appreciate the benefits of partial CSI at the transmitter in the proposed repetition protocols, 
consider the performance with perfect CSI (F = +oo). 



Since the power control is causal, the throughput is found as a solution of a dynamic program [13|: 



Proposition 9. The throughput vi^m f=oo ib '■^ ^^^^ solution of: 



1 - PrE£i Ut < gm 



max R- 



s.t. 



E 



M 

m=l 



E 



P 1 



{J:7Li Ut<g{R)} 



1 + E':!:'! pr[Er=i Ut < 9{R)] 



(39a) 



(39b) 
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-{7m-Pm>e«-l} 



indicates successful decoding on the m-th transmission and g{R) 



where: for ALO, Um 

0;forRTD, Um = ImPm represents the received SHR on the m-th transmission and g{R) = > for 

INR, Um = log(l +7m-Pm) represents the mutual information at the receiver on the m-th transmission and 
g{R) = R; the optimization is with respect to causal power policies Pm = Pm{li, ■ ■ ■ , 7m), fn = 1, . . . , M. 



t-i 



Proof: The optimization of a causal pow er control system can be cast as a dynamic program over 
finite horizon with complete observations [ [l3| , |15|, where the system state is = ^ 
control is Ut, and the state evolves as Xt+i = Xt -\- Ut, with Uq = 0, from which (|39l) follows. 



Um^ 



the 



The problem in P9a[ ) is similar to the outage minimization problem in [15|. As in [15|, we can write 
the iterative algorithm that defines the optimal dynamic programming solution, however an explicit closed 
form solution is not available in general. Numerical techniques, as those proposed in p3| , must be used 
for numerical evaluations of (iBQl). 



VIL Novel bounding technique 

In the previous section we proposed novel protocols that combine power control with partial CSI and 
retransmissions. In all cases, a closed form expression of the throughput requires a closed form expression 



for the cumulative density function of: (a) the fading 7^ for pl^Y^^ in ([36]), (b) random variables of the 



011 



type Sm = J2T=i 7^' f*^^ I^m^^^ ( |37] ), and (c) random variables of the type = 11^=1 (^^ 

for some > for p^^^^ in ([38]). Since the distribution of is rarely known in closed form, in the 
following we propose a novel bounding technique for the cumulative density function of in terms 
of the the cumulative density function of Sm- For the case of iid Rayleigh fading, the density of Sm is 
known in closed form; hence, our technique allows to determine closed-form upper and lower bounds for 
probabilities involving /„. 

Consider generic non-negative constants {Tm}, m G N, a constant > 0, and define 



„(ALO) 
fm 



Pr 



7. 



< 1, s 



l,...,m 



.s=l ' 

— ^Vlogfl + e^) <1 

for some sequence {7^} of iid random variables. 
Example: As an example, consider Fig. [sj which shows in the plane Xi 



(40a) 
(40b) 

(40c) 



defines -k G {ALO, RTD, INR}. The probability ^2^^"^'' is the integral of the joint density of (xi, X2) 
over the square {xi,X2) G [0, 1]^. The probability p'^^^^ is the integral over the triangle xi + X2 < 1, 
a;i > and X2 > 0, that is, over the region in the positive quadrant below the dotted-line curve labeled 
"RTD" in Fig. 5 And finally, the probability P2^^^ is the integral over the region in the positive quadrant 
below the solid-line curve labeled "INR" in Fig. ]5] The curve labeled "INR" in Fig. ]5]is a convex function 
that can be bounded from above and from below by piece-wise linear functions. We chose piece-wise 
linear functions because the region they define is the union of triangular regions. In particular, for the 
inner bound, we take the union of the two regions below the tangent lines at (xi,X2) = (0,1) and at 
(xi, X2) = (1, 0) (the region in the positive quadrant below the dash-dotted-line curve labeled "INR inner 
region" in Fig. ]5]), while for the outer bound, we take the union of the two regions below the lines passing 
through {xi,X2) = (0, 1) and {xi,X2) = (1/2, 1/2), and through (xi, X2) = (1,0) and (xi, X2) = (1/2, 1/2) 
(the region in the positive quadrant below the dashed-line curve labeled "INR outer region" in Fig. ]5]). 
By extending the idea presented in the above example to the case of a general m G N we can show: 



(ALO) 



^ , X2 = ^ the regions that 
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Fig. 5. The region that defines the probability of outage for INR with M — 2 and its approximations. 



Proposition 10. The probability in (40c) can be bounded as: 



(INR) y p 

^'^^^ < Pr 



7, 



7i 



s=l 



9 max — < 1 

Tg t=l,...,m Tt 



(41a) 



; — n m max — < 1 

;i + ey/"' - 1 J t=i,...,m ^ 



(41b) 



Proof: The proof can be found in AppendixjD ■ 
The interesting fact about the two bounds in is that they are computable from the knowledge of 
the density of the random variable 



Xa b = CL max — + b 

t=l,...,m Tt 

s=l 



(42) 



for some fixed (a, 6) G M^. 

Proposition 11. For the case ofiid negative exponential random variables {7t} f/.e., iid Rayleigh fading), 
the density of Xa^b in ( |42| ), for any (a, 6) G M^, /5 g?ve/i m (44) in Appendix^ 



Proof: The proof can be found in Appendix lEl 



VIII. The iid Rayleigh fading channel 



To illustrate the gain achievable with the protocols proposed in Section VI we evaluate the performance 
of the different protocols for the Gaussian iid Rayleigh fading channel, for which the fading cumulative 
distribution function is Fy(x) = 1 — e^"^' for x > 0. We define the exponential integral function as: 



E 



7 



RF 



POO 

/ e-'/tdt^Ei{x) 

J X 



for a; > 0. We use the symbol 



iiRF,, 



to indicate that the equality holds for Rayleigh fading channels. 
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Fig. 6. Ratio between tlie throughput and Va^J^ f=oo ('•'^^ ergodic capacity with full CSI), for the Rayleigh fading channel. 
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Fig. 7. Ratio between the throughput and Vm^^ f=oo (^^^ ergodic capacity with full CSI), for the Rayleigh fading channel. 



For the plots, communication rates are measured in bits/sec/Hz and the figures show the relative 
throughput performance with respect to the ergodic water-filling capacity (i.e., M = F = oo and INR), 
which is the ultimate performance limit for a fading channels with full CSI. Table |l] summarizes the cases 
considered in the following. 



TABLE I 



Amount 


M 


= 1 


M 


= 2 




M 


= oo 


of CSI 


outage cap. 


HARQ protocols 


erj 


;odic cap. 


Absent 


in 


VIII-Bl 




(impossible, need at 


in 


VIII- A 1 


(F = 1) 






least 1-bit for ack/nack) 






1 bit 


in 


VIII-B2 




in 


VIII-C 


(classical HARQ) 


in 


VIII-A2 




(F = 2) 






in 


Vlll-D 


(proposed HARQ) 






Full 


in 


VIII-B3 




(not evaluated) 


in 


VIII-A3 




(F = oo) 

















A. Throughput upper bound (ergodic capacity) 
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(INR) RF i„. 



1) M = oo, F = 1, INR: With constant power allocation the ergodic capacity is: 

1 

2) M = oo, finite F, INR: The ergodic capacity with partial CSI is given in Theorem |4| 

„(INR) 
RF 



max { 5^ - log(l + PfSf)e-'f - e^/^fEi{l/Pf + Sf) 
+ log(l + PfSf^i)e-'f-^ + e^/^fEi{l/Pf + 



where the thresholds {s/} are defined in (|7]), the maximization over {P/} is subject to the constraint 
in ^ and 



F-l 

RF 



The ergodic capacity with partial CSI can be bounded by using Proposition |5| in particular, Hf in ( |13a| ) 
is given by 



RF e-^/-i(l + s/_i) -e-^/(l + s 



f) 



3) M = oo, F = oo, INR: With water-filling power allocation the ergodic capacity is: 

^E,F=oo = Ei(A), 

pM^-Ei(A). 
A 

4) Discussion: From Fig. [g] we see that rj^^^^ (i-^-' ^ bit of feedback) is already at 97% of the 
water-filling capacity ?7^^^^^^^ at an SNR as low as -25dB. At high SNR, rifj^^2,,F=i i^*^ CSI) behaves 
like rjfj^^^p^^ (full CSI) but at low SNR their ratio tends to zero. This restates a well know fact that 
power allocation in single user channels offers benefits at low SNR only and that a single bit of feedback 
(i.e., F = 2) gives almost all the gain achievable by full CSI (i.e., F = oo). 

B. Throughput lower bound (outage capacity) 

1) M = 1, F = 1, ALO: With constant power allocation the outage capacity is: 

VMt%=i ^ maxlog(l + Psi)e-'\ 

' si>0 

2) M = 1, F = 2, ALO: With 1-bit of feedback, the throughput is the solution of: 

(ALO) _ -si 

Vm=i,f=2 — „ , max e 

0<si<.S2<oo 

log I 1 ^ 



S2 Si S2 

max e"'*' log(l + Psie'"^). 

0<si<oo 
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3) M = 1, F = oo, ALO: With truncated channel inversion power allocation the outage capacity is: 



(ALO) RF 1 A , ^ ^ 

= ina^log + 



Vm 

4) Discussion: When M = 1, there is no need to send a ACK/NACK because the transmitter cannot 
retransmit. In this case, the one bit of feedback should indeed be used at the beginning of the slot to 
inform the transmitter about the state of the channel. From Fig. |6]we observe that 1-bit of feedback at 
-25dB results in ri[^t%=2/vi^SL,F=oo = §1% while (with constant power) ?7jt=?F=i/^Ki,F=oo = 10% 

only. At +5dB, ri[^t%=2/VAl^=l,F=oo = 73% while (with constant power) ri[^ti^p^i/VA^^=l,F=oc = 45%. 
In fact, at high SNR, power allocation is less critical and the gain due to CSI vs. no CSI diminishes. 
We also reported for comparison the achievable throughput for ALO with F = 4 (2 bits of feedback) by 
using the approximation S£ = si^^~^ as in Proposition [v} and for ALO with F = oo (full CSI). We see 
that the gains attainable at low SNR due to only a few bits of feedback are dramatic and that 2 bits of 
feedback attain a throughput remarkably close to the case with full CSI. 



C. Classical HARQ protocols with F = 2 and M = 2 

Classical HARQ protocols use the 1 bit of feedback (i.e., F 



2) to signal ACK/NACK. Classical 



HARQ protocols are a special case of the protocols proposed in Section VI obtained by setting Smj = oo 
for all m > and / > in ( [25] ). In classical HARQ protocols the power can vary across repetitions (i.e., 
different values of Tm, m = 1, . . . ,M) but it cannot depend on the CSI [7|. In this case the throughput 
in ( [28] ) is a function only of {pm,o}, me {1, . . . , M + 1} given by: 

m 

11(1-^'") 

1 - e-^= 



~(ALO) RF 
Pm,0 



s=l 
m 



~(RTD) RF 
Pm,0 ~ 



J^INR) 



Pr 



Y^\og[i + e^-^\ <iog(i + 

s=l ^ ^ 



for 6 = 



1 > and m > 1. Consider the case M = 2: the only probabilitv not known in c losed form 



is Pm^(^\ which we bound by using the technique developed in Proposition [ll| in Section 

1 



VII 



as follows: 



—(in) 
^2,0 



< ^INR) 
P2,0 



^ ~{out) 
< P2,0 



^ ~(RTD) 
<F2,0 

^ /'2,0 



1-g 



E 



9 



7iM 



2 + e 
1 

""'1 + 71/^0;] 



1-9(1/2) = 
(l-e--)(l 



I-T-i/t-2 \-T2lTx 



where the function q{x) is given by (|44]) in Appendix |e| 
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1) Discussion: In Fig. |7| classical HARQ protocols are labeled as "classical" in order to distinguish 
them from the novel protocols proposed in this work, which are labeled as "new". The case of M = 2 
transmission is considered (which implies that the throughput needs to be optimized with respect to the 
two parameters (r2,ri)). From the numerical results for SNR = P E [— 25, +25]dB, we saw that equal 
power allocation across transmissions (i.e., T2 = ri) is optimal at all SNR's for ALO; for RTD and INR, the 
use of different power (i.e., T2 ^ n) offers benefits; however the improvement is negligible (for example 
less than 0.3% across the entire range of simulated powers for RTD) for this reason in Fig. |7] we only 
show the throughput with the optimized {t2,ti) only for INR. We see that the throughput trend for the 
classical HARQ protocols with F = 2 and M = 2 is the same as the one for F = 1 and M = 1 reported 
in Fig. [6j This shows that the 1 bit of feedback used for ACK/NACK only does not offer substantial 
throughput improvement compared to the outage capacity with constant power allocation (F = 1 and 
M = 1). Only for SNR's larger than lOdB, using classical INR with F = 2 and M = 2 gives a larger 
throughput that ALO with F = 1 and M = 1 with power control. From these observations we conclude 
that classical HARQ make an inefficient use of the feedback resources. 



D. New protocols for M = 2 and F = 2 

Since we consider the case F = 2, we only need to characterize Pmj, for m G {1, 2, 3} and / = 0, as 
per Proposition [8} 

The probability of requesting a retransmission is: 



Pi,o = Pr[5i = 0] = Pr 



7i 

Sl,l 



< 1 



Pr[r = 2] 



the probability of decoding failure after the first transmission is: 



P2,0 


= Pr[fii 


= 0, 


B2 = 


= 0] = 


Pr 


71 




<1, 


■52,1 


min{ri, 


Sl,l} 


Pr 


71 




<1, 


7i ^ 72 

Tl S2,l 


min{ri, 


ri,i} 


Pr 


7i 




<1, 


7i ^ 72 


min{ri, 


n,i} 





< 1 



'-(i + e 
1 ^ 



< 1 



ALO 
RTD 
INR, 



and the probability of decoding failure after the second transmission, which coincides with the probability 
of outage, is: 



Pr 
Pr 
Pr 



= Pr[5i 


= 0, 


7i 




min{ri, 


■51, i} 


7i 




min{ri, 




71 





<1, 



<1.^+ 



72 



<1 



ALO 



72 



<1 



RTD 



72 



1 + 9 



71 



<1 



INR. 



_min{ri, ' Si min{r2, S2,i} v n. 

The probabilities for ALO and RTD can be evaluated in closed form, while the probabilities for INR 



can be bounded as in Subsection |VIII-C (from Proposition [TT] in Section |VII[ ) for the classical HARQ 
protocols by using the function q{x) in (|44l). 
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1) Discussion: Fig. |7] shows the ratio among the throughput of new protocols and the ergodic water- 
filling capacity. We see that the new protocols dramatically outperform the classical repetition protocols, 
especially at low SNR (compare with Fig. [6]). Indeed, at low SNR it is critical to be able to save power 
when the channel is in deep fade. By providing the transmitter with a 1-bit quantization of the current 
channel gain (rather than ACK/NACK), we enable the transmitter to do so. At low SNR, the repetition is 
not needed (ALO with M = 1 and F = 2 has the same throughput of INR with M = 2 and F = 2), while 
at high SNR, the repetition helps. Notice that here high SNR means is SNR > 5dB, at which 'ini^^^2 f=2 
about 67% of r^J™^^^^. We did not report the RTD and ALO curve for M = 2 and F = 2 with CSI 
as they do not differ much from ALO with M = 1 and F = 2. 



IX. Conclusions and Future Work 

In this work we considered HARQ protocols where the feedback bits not only convey a retransmission 
request to the transmitter but also inform the transmitter coarsely about the channel state. We developed a 
new class of protocols that feedback the quantized index of a suitably scaled version of the current fading 
value; the scaling factor is such that the mutual information already accumulated at the transmitter from 
the previous transmissions is taken into account. We showed that our proposed protocols significantly 
outperform classical HARQ protocols for the same amount of feedback resources, especially at low SNR; 
this shows that ACK/NACK feedback is suboptimal in time-varying channels. 

As future work, it would be interesting to evaluate the throughput performance when the cost of 
acquiring the CSI and the error in the estimated CSI are taken into account. Also, it is important to test 
the proposed protocols with practical codes, instead of with ideal Gaussian codes. 



Extensions to multiple access channels are presented in [18|. 
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Appendix A 
Proof of Theorem [T] 

As a direct application of the renewal-reward theorem f2)\, we have that the limit in (|3]) converges almost 
surely to E[P]/E[7~] and the limit in Q converges almost surely to E[7^]/E[7~]. The average inter-renewal 
time E[7^ is given by: 

M M 

E[T] = Pr[r = m] = Y^ Ft[T > m]. 

m=l m=l 

The average reward E[7^] is given by: 



E[7^] = R{1 - Pout) = R{1 - Pr[r = M, fail to decode]). 
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since in our framework, a packet is lost (TZ = 0) only if on the last transmission decoding was not 
successful. The average transmit power E[V] is given by: 



AI 



E[V] = Pr[r = m]E[V\ T = m] 

m=l 

M m 

= ^Pr[r = m]^E[Pt|r = m 

m=l t=l 

M M 



t=l m=t 
M 



E^=tPAr = m]E[P,\r = m] 



M 

Y^PT[r>t]E[Pt\r>t], 



t=i 



QED. 



Appendix B 
Proof of Theorem [6] 

Given a partition {TZf} of IR+, the transmit power is defined through the thresholds in (|7]) as 

/=i ^ 

and must satisfy the average power constraint 

A e« - 1 - 



'7 



For a fixed rate i?, the throughput is maximized if the outage probability is minimized. The outage 
probability satisfies 



Pr 



e«- 1 



{7e7^/} 



< P 



/=i 



E 



E 



{log(l+7Ef=l^l{.67^,})<R} 



l{7e7e/}l{7<s/} 

./=i 



> E 



min 1 



where the last inequality holds with equality for 



{7<s/} 



= {7 e : / = arg^min^ 
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By assuming without loss of generality that the thresholds {s/} are ordered in increasing order, we have: 



are; min 1 



{7<s/} 



G {1,...,/} if 7 e / = 1, 

G {!,..., F} if 7 e [SF, +oo) U [0,si) 

In order to use the least power we choose 

are min 1; 



-{7<s/} 
,(ALO) 



= / if 7 e i?f = [S/, S^+l), / = 1, . . . ., F - 1, 

= F if 7 G R^'-'^^ = [sF, +oo) U [0, si). 

With these quantization regions, an outage only occurs when the feedback value is F and the fading is 
in [0,si), hence the throughput is 

„(ALO) 



max R Pr 



max log I 1 + 

{Sf} 

■ Pr[7 > si] 



i?< log I 1 + 7X1^^1 



P 



Appendix C 
Proof of Proposition [7] 

^ALO) 



The optimal values ofO<S2<...<Si7'< sf+i = oo for rj)^^i f in p2| ) minimize 



[^7(^/+l)-^7(^/)]- 



(43) 



By assuming that the fading has a density f^{x) = dFj{x)/dx, by taking the partial derivatives of ( |43| ) 
with respect to si for £ > 2, and solving them equal to zero, we get that the thresholds satisfy 

se-i Sify{si) ' ~ 

For sufficiently large F, since the thresholds are going to be close to each other, we can approximate 

F^(se+i) - F^{se) dF^(x) 



(s^+i - Si) dx 
and hence we conclude that the optimal thresholds satisfy 

s'j ^ S£_iS£+i Si ^ Si ^^"^ 

for some ^ > 1. 
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Appendix D 
Proof of Proposition [TOl 

For a general m G N, the probability in ( |40c[ ) (with = is equivalent to 

Pr[X™-/(Xi,...,X„_i) <0], 



„(INR) 



for 



/(xi, . . . ,a;m_i) = \ 



l + d 



Our goal is to bound the region below /(xi, . . . , x^-i) in the positive quadrant by regions defined as union 
of hyperplanes. Toward this goal, it is useful to keep in mind that the partial derivatives of /(xi, . . . , Xm-\) 
are given by 

df{Xx, Xm-l) l + O /(Xi, . . . , Xra-l) 



dxi 



1 + 0X, 



Being /(xi, . . . , x^-i) a concave function in all its arguments, it can be lower-bounded by any hyperplane 
tangent to it. In particular, if we consider the hyperplanes tangent to /(xi, . . . , Xm-i) at those points with 
at most one non-zero coordinate at value 1, we obtain: 



/(xi,. . . ,x„_i) > 



for all t = 1, . . . , m, which is equivalent to 



= l,...,m— 1, s^t 



/(xi,...,x^_i) < 



l + xtO 

i + e ' 



, m. 



This bound implies 



(INR) > p 



U(i+<')E7<i + 

t=l s=l 



Pr 



1 + 0) V — < 1 + ^ max 

^ r, t=l,...,r 



7t 

m Tt 



( |4Ti| ). 



In the same spirit, we bound /(xi, . . . , x^-i) from above by considering the union of the region below 
the m hyper-planes defined as follows: for a fixed t G {1, . . . .m}, the hyper-plane that passes through 
the points with coordinates x^ = 1 and Xj = for all j ^ k and 7^ t (this defines m — 1 points with 

^ — '-7, IS defined as 



only a single non-zero coordinate at 1) and the point xi 



The region below the union of the above hyperplanes for alH = 1 
p^^^ and hence 



m ] Xt = 1- 

, . , m contains the region that defines 



'mm , Q V 



Pr 



m I max — < 1 

t=l,...,m 



(|4Tbl). 
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Appendix E 
Proof of Proposition [TT] 

It is a well known result in order statistics (W, Ch.5] that the order statistics of a sample of independent 
negative exponential random variables can be expressed as the unorder statistics of a sample of independent 
negative exponential random variables with appropriate mean value. We report the derivation of result 
here for sake of completeness. 

Let {7fc}fcLi be an independent sample of size K of negative exponential distributed random variables 
with mean E[7fe] = Hk- Let tt be a permutation of the integers {1,2,..., K] and let Vk be the set of the 
K\ such permutations. With an abuse of notation, we also indicate with tt the event 

A 



It is well known that 



Pr[7r] = Pr[7^(i) > 7^(2) > . . . > 77r(A-)], Vtt G Vk- 



1^1 X 

;m n ~ ( ~ ~ ) i{-.(i)>->-.(x)>o}- 



Fix a permutation vr, define '^■^{k+i) = 0, and consider the following change of variables with c-'^^ > 
for alH = 1, . . . , i^: 

K 



Zi = cf ^(7^(i) - 7^(*+i)) Mi) = 5Z 



The random variables {ZkY^^i non-negative by definition. Moreover, the transformation giving {Zk]^=i 
from {7fc}|Li has Jacobian HiliCi''''- Hence, the joint density of {Zk}^^^ is 



fzi,...,ZK\n{Zi, . . . , Zn) 

K K 



1 1 

exp 



Prfvrl „. 
^ 1 




r- 



tier' ^ 9) 



where 



(tt) 



obtained by recalling that 



K / i 



Pr TT . ^ 



HE 



/ivrO-) 



^ /^7r(fc) 



In words, the change of variables has produced the unorder statistics of an independent sample of size K 
of negative exponential distributed random variables with mean E[Zfe] = e^^\ Remark: in the case fii = fi 
for alii = 1, . . . , one obtains the familiar result (think at the inter- arrival times of a Poisson process): 

>) j' Pr[7r] 
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We now apply this result to the computation of the distribution of 



K 



k=l 



Ik 



For general a and h we have: 



Pr 



K 



a max^7fc + & 7fc < 

k=l 



X 



^ Pr[7r] Pr 



K 



a77r(l) + ^ ^l^{k) <X\'K 



k=l 



a; TT 



^=1 



= ^ Pr[7r] Pr 



Elr + 'E^^s 



X TT 



We next chose 4'''' = ^^^^ so that the are iid negative exponential with unit mean value for 
all ^ e {1, . . . , K} and for all tt e T'x, hence we obtain 



1 



€=1 



e=i 2^k=i ij.^. 



(k) 



At this point, the problem reduces to that of finding the density of a random variable Y that is a linear 

combination of iid negative exponential with unit mean random variables, i.e., Y = VkZ^ for Vk 

for all /c = 1, . . . , . The calculus of residues applied to the characteristic function of Y, assuming that 
all the coefficients Vi are distinct, gives: 

K 



k=l 



K 



= JjE[exp(-ja;vfcZfe)] 

k=l 

K ^ K 

— TT — \^ Q'fc 

with ttfe = Residue [fl'y (a;), = ]^ ^ _ , , 

K 

such that ^ afc = 1. 



K 



k=l 



By taking the inverse Fourier transform of gvi^^) we obtain the density function of Y given by 



frix) = ri"^ "^^"^ l{xsign(^,)>0} 



it=l 
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and hence 



k=l 



fe=l;'i;fe>0 
K 

+ [ J2 «,e-^/"'=] l{x<o} 



k=l:vk<0 



If some coefficients are equal, say vi and V2, then it suffices to late the limit for vi — > V2 of Pr [ X^i^i ^fe-^fc < 

Back to our original problem: 

Pr[X„,,<a;]= ^ Pr[7r][(l- J] «f V^/"^"') l{,>o}+ 



M -x/v 



(n) 



ai e 



4a;<0} 



Vi = 



(tt) a 



^ 1 

ii ^1 (tt) / (tt) 



For example, with 6 > and a + 6 > 0, so that a + bi>0 for all £ e N (this is the case of interest in 
our problem), for K — 1: Xa,b — {a + b)^i, for x > 

For K ^2: Xa,h = (a + 6) max{7i, 72} + 6min{7i, 72}, for x > 



-XTl 



CT2 



+ e" 



-XT2 



CTi 



_j_ q-Cx{ti+T2) 



cri + CT2-ri CT1 + CT2-T2 

TiT2(l - 2c) 



+ 



1 



{CT1 + CT2- ri)(cTi + CT2 - T2) ' 

a + b 



[a + 6)/Xjk a + 26 



(44) 
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